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Orbital fluctuation theory in iron-based superconductors: s++-wave superconductivity, 
structure transition, and impurity-induced nematic order 
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The main features in iron-based superconductors would be (i) the orthorhombic transition ac- 
companied by remarkable softening of shear modulus, (ii) high-Tc superconductivity close to the 
orthorhombic phase, and (iii) nematic transition in the tetragonal phase. In this paper, we present 
a unified explanation for them, based on the orbital fluctuation theory, considering both the e-ph and 
the Coulomb interaction. It is found that a small e-phonon coupling constant (A ~ 0.2) is enough 
to produce large orbital (=charge quadrupole O^z/yz ) fluctuations, which causes the s-wave su- 
perconductivity without sign reversal (s++-wave state). The derived orbital fluctuations also cause 
the instability toward the structure transition due to the bound state formation of two orbitons 
with opposite momenta, which is called the "two-orbiton process". Moreover, impurity-induced 
non-local orbital order with C2-symmetry is obtained when the orbital fluctuations are strong. This 
"impurity-induced nematic state" explains the in-plane anisotropy of resistivity in detwinned sam- 
ples. We stress that (i)-(iii) are reproducible only when orbital fluctuations with respect to Oxz 
and Oyz charge quadrupoles are the most divergent. This fact ensures the reliability of the present 
model Hamiltonian and calculation. 
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I. INTRODUCTION 

Iron-based high-Tc superconductors had been discov- 
ered by Kamiahara et al. in 2008^, and the highest su- 
perconducting (SC) transition temperature T^ at present 
reaches 56K. The SC state occurs when the crystal struc- 
ture is tetragonal, which is realized by chemical dop- 
ing (or applying the pressure) on mother compounds, 
through the structure transition from orthorhombic to 
tetragonal. In the phase diagram, the SC phase is next 
to the non-SC and metallic orthorhombic phase, and 
the structure transition a,i T ~ Ts is second-order in 
Ba(Fe,Co)2As2^. Very large softening of shear modulus 
Cs suggests the existence of strong orbital (quadrupole) 
fluctuations^"^. Moreover, the spin-density- wave (SDW) 
state with Q sa (tt, 0) occurs in the orthorhombic 
phase, that is, Tn is close to but always lower than 
Ts- These experimental facts suggest a close relation 
between the mechanism of superconductivity and struc- 
ture/orbital/SDW transition. 

Although the lattice deformation in the orthorhombic 
phase is very small ((a — h)/{a + h) < 0.003), in-plane 
resistivity shows sizable anisotropy in the orthorhom- 
bic phase. This fact means that the structure transi- 
tion is driven by electron-electron (or electron-optical- 
phonon) interaction, not by the cooperative Jahn- Teller 
effect due to electron-acoustic-phonon interaction. Very 
interestingly, large in-plane anisotropy starts to occur in 
the tetragonal phase at T*, which is about lOK^lOOK 
higher than Ts in detwinned Ba(Fe,Co)2As2^. The dis- 
covery of this "nematic electronic state" free from lat- 



tice deformation had attracted great attention. The ne- 
matic transition also occurs in the tetragonal phase of 
BaFe2(Asi_2;Pj,)2, confirmed by the in-plane anisotropy 
in the magnetization (xa 7^ Xb) using the torque measure- 
ment under magnetic field^. These experimental facts of- 
fer us great hints to understand the electronic states and 
the pairing mechanism in iron-based superconductors. 

Now, we have to try to construct a theory that can 
explain abovementioned main characters in pnictides in 
a unified way, not restricted to the superconductivity. In 
the early stage of the study of iron-based superconduc- 
tors, however, many theorists had concentrated on the 
study of pairing mechanism. Based on spin fluctuation 
theories, fully-gapped sign-reversing s-wave (s±-wave) 
state had been predicted^"^^. The origin of the spin 
fluctuations is the intra-orbital nesting and the Coulomb 
interaction. However, the robustness of Tc against ran- 
domness in iron pnictides indicates the absence of sign- 
reversal in the superconducting (SC) gap^'^"-^^. Later, 
orbital-fiuctuation-mediated s-wave state without sign 
reversal (s++-wave) had been proposed^^"^^. The orbital 
fluctuations mainly originate from the inter-orbital nest- 
ing and the electron-phonon (e-ph) interactions due to 
non-Aig optical phonons. 

One of the main merits of the orbital fluctu- 
ation scenario is the robustness of the s++-wave 
state against impurities or randomness, consistently 
many experiments^"*"^^. Moreover, orbital- fluctuation- 
mediated s++-wave state scenario is consistent with the 
large SC gap on the 2^-orbital band in Bal22 systems^®, 
observed by bulk-sensitive laser ARPES measurement^*^. 
Note that the "resonance-like" hump structure in the 
neutron inelastic scattering^ is frequently explained as 
the spin-resonance due to the sign reversal in the SC 
gap^^'^^. However, experimental hump structure is well 
reproduced in terms of the s++-wave SC state, rather 



than the s±-wave SC state, by taking the suppression in 
the inelastic scattering 7(0;) for \oj\ < 3A in the SC state 
(dissipationless mechanism)^^'^"'. 

In this paper, we present recent developments of the 
orbital fluctuation theory presented in Refs.^^"^^'^^'^®. 
We present a unified explanation for the the following 
main characters in iron-based superconductors: (i) or- 
thorhombic transition accompanied by remarkable Cs 
softening"^"^ , and (ii) emergence of high-T^ superconduc- 
tivity strong against randomness next to the orthorhom- 
bic phase. We also discuss (iii) "nematic transition" in 
the tetragonal phase^ in terms of the impurity-induced 
non-local orbital order. It is noteworthy that (i)-(iii) can 
be explained only when orbital fluctuations with respect 
to Oxz and Oyz charge quadrupoles are the most diver- 
gent. This fact assures the reliability of the present the- 
ory and model Hamiltonian. 



II. ORBITAL FLUCTUATIONS AND S++-WAVE 
SUPERCONDUCTIVITY 



A. Antiferro-orbital fluctuations due to 
inter-orbital nesting 

Here, we study the five-orbital tight-binding model in- 
troduced in Refs.*'^ , which reproduces the experimental 
multiband structure very well. We also include both the 
Coulomb interaction ([/, U' , and J) and the quadrupole- 
quadrupole interaction induced by optical phonons. In 
this paper, we introduce the xy2:-coordinate, in which x- 
and j/-axes are along Fe-Fe direction. 

Because of the symmetry of AS4 tetrahedron, Fe- 
lon optical phonon induces the following quadrupole- 
quadrupole interaction:^*, 



(r 



When g = 0, five quadrupole susceptibilities Xp(q) 



xz,yz,xy, z^,x^ - 



y ) induced by Coulomb interac- 



tion are small and almost comparable, with peaks at q = 
(tt, 0) or (0,0). By introducing small quadrupole inter- 
action g, quadrupole susceptibilities with F — xz,yz,xy 
are largely enhanced. The most divergent susceptibilities 
are X^zCq) at q = (tt, 0) and x%{^) at q = (0,7r): They 
are approximately expressed as 



X^q) 



ce 



l+e2(Qr-q)2-2Cc>/c^o 



(2) 



for T — xz and yz, where Q^z — {t^, 0) and C^^z — (0,pz). 
In two-dimensional systems, orbital fluctuation parame- 
ters behave as^^ 



e - i{T-TAFQr 

luq = l'(T - Tafq) 



(3) 
(4) 



The Weiss temperature Tafq becomes zero at orbital 
quantum-critical-point (QCP). In eq. (3), the unit of 
the length is aFe-Fe(~ SAngstrom). 

In later sections, we will explain that the develop- 
ment of xiz(yz)^^xz(vz)) gi'^^s rise to not only the s++- 
wave state, but also the structure transition as well as 
impurity-induced nematic transition in the tetragonal 
phase. Although superconductivity can be caused by 
other orbital fluctuations (such as F = xy, x^ — y"^, or 
z^), both structure and nematic transitions are uniquely 
explained only when T = xz and yz are the most diver- 
gent. Therefore, the present orbital fluctuation model is 
considered to be reliable. 



B. Phase diagram for s±- and s++-wave states 
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where Or,i is the quadrupole operator for channel F at 
site i introduced in (I). g{uJi) = g{0) ■ uj'^/{ujf + cj|j) is 
the quadrupole coupling induced by the optical phonons, 
were ujjj ~ 200 ~ 300K is the phonon frequency and 
Wi = I-kTI is the boson Matsubara frequency. Note 
that Opi/ ex /^/^ + /^/^, where I is the angular mo- 
mentum. Here, we set (^xz\Oyz\xy) — {yz\Oxy\xz) ~ 
{xy\Oxz\yz) = 1 by multiplying a constant. Recently, 
we have found that Eq. (1) is also caused by the mul- 
tiorbital Coulomb interaction, by including the multior- 
biton exchange process' that is absent in the random- 
phase approximation (RPA) [S. Onari and H. Kontani, 
arXiv:1203. 28741]. 

In iron pnictides, antifcrro-quadrupole (AFQ) fluc- 
tuations with respect to F = xz/yz are induced by 
the quadrupole interaction and the inter-orbital nest- 
ing. The quadrupole susceptibility is Xp(q, t) — 

/o''dTe^"'-(r,6r(q,T)Or(q,0)). 



Figure 1 (a) shows the phase diagram obtained by the 
RPA (mean-field approximation) for the electron filling 
n — 6.1, where 17 (0) is the quadrupole interaction and 
U is the intra-orbital Coulomb interaction; we set the 
Hund's coupling term J — C//10 and the inter-orbital 
term U' — U — 2J. Uc (««) is the charge (spin) Stoncr 
factor; ac — I {as ~ 1) corresponds to the orbital (spin) 
ordered state. For U — 1.3eV, orbital order appears 
when g{0) « O.lGeV, meaning that about 50% of orbital 
fluctuations originate from the Coulomb interaction. By 
solving the SC gap equation, we obtain the s++-wave 
state when the orbital fluctuations dominate the spin 
fluctuations^^. 

In Fig. 1 (b), we show the U-g phase diagram ob- 
tained by the FLEX approximation. The dashcd-dotted 
line represents the condition ac — 0.98 at T = 0.015, 
corresponding to 17 = 0.25 ^ 0.3. Therefore, substantial 
orbital fluctuations emerge for A — gN{0) < 0.2 even if 
the self-energy correction is taken into account. On the 
other hand, a^ == 0.95 (0.92) for t/ = 1.8 and 5 = 
(0.3) in the FLEX approximation, although C/cr — 1-25 
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FIG. 1: g-U phase diagram obtained by (a) RPA for J /U = 
1/10 and (b) FLEX approximation for J/U = 1/6. Also, 
(c)-(e) show the SC gap structure in FSl (inner h-pocket), 
FS2 (outer h-pocket) and FS4 (e-pocket) for different model 
parameters. Sk ~ tan~^((fcy — ky)/{kx — k^)), where k° is the 
center of each FS. 



for as = 1 in the RPA. Thus, the renormalization in a^ is 
rather larger than that in ac- The region of the s++-wave 
state is widely expanded in the presence of dilute impu- 
rities. In the present model, both x^^i z) ^^'^ ^xy ^^^ 
strongly developed. Since all t2g-orbitals on the Fermi 
surfaces are involved in these fluctuations, we obtained a 
fully-gapped s+-|--wave state as shown in Fig. 1 (c), con- 



sistently with many high-Tc pnictides. When the spin 
fluctuations are comparable to orbital fluctuations, on 
the other hand, competition of these fluctuations induce 
the nodal structure in the SC gap, as shown in (e). 

In summary, we discussed both the spin-fluctuation 
mediated s±-wave state and orbital-fluctuation mediated 
s+4.-wave state in iron-based superconductors. We have 
shown that the latter is realized by inter-orbital nest- 
ing in the presence of small quadrupole interaction g in 
eq. (1). Figure 1 (b) show the U-g phase diagram ob- 
tained by the FLEX approximation. Since both s±- and 
s+-|--states belong to the same Aig symmetry, a smooth 
crossover from s±-state to s++-state is realized as in- 
creasing g or impurity concentration riimp. During the 
crossover, nodal s-wave state is realized, at which the 
obtained Tc is rather suppressed but finite. (Mathemat- 
ically, it is impossible to distinguish between nodal s++- 
and nodal s±-states.) Experimentally, nodal line struc- 
ture appears in Ba(Fei_a;Coa;)2As2, in which both spin 
and orbital fluctuations are considered to be developed^°. 
Recently, we have verified that the realization condition 
for the nodal structure is much moderate if we study a 
realistic "three dimensional" modeP^. 



C. heavily hole-doped system: KFe2As2 

Here, we study the superconducting state in newly dis- 
covered heavily h-doped superconductor Ka;Fe2Se2 {Tc ^ 
30K) based on the ten-orbital modeP^. Since the hole- 
pockets are absent, the s-t-wave state is unlikely. How- 
ever, when the Coulomb interaction is large enough, spin- 
fluctuation mediated d-wave state would appear due to 
the nesting between electron-pockets^^. However, the 
symmetry of the body-centered tetragonal structure in 
K2:Fe2Se2 requires the existence of nodes in the d-wave 
gap^*^, although fully- gapped d-wave state is realized in 
the case of simple tetragonal structure. In the presence of 
moderate quadrupole interaction g, on the other hand, we 
find that orbital fiuctuations give rise to the fully-gapped 
s++-wave state. 

In Fig. 2 (a) , we show the ac-dependence of the eigen- 
value of the gap equation, Ae, for the s++-wave state 
with [7 = 0, and the ag-dependence of Ae for the d-wave 
state with g =■ 0. In calculating the s++-wave state, we 
use rather larger phonon energy; wd = 0.15 eV, consid- 
ering that the calculating temperature is about ten times 
larger than the real Tc. The SC gap functions for s++- 
wave state are shown in Fig. 2 (b). The SC gap would 
become more anisotropic in case of f7 > 0. 

Therefore, similar to iron-pnictide superconductors, 
orbital-fluctuation-mediated s++-wave state is realized 
by small e-ph interaction A = gN{0) ~ 0.2. The obtained 
Ae for the s_|_+-wave state is larger than that for the spin- 
fiuctuation-mediated d-wave state. We stress that nodal 
lines appear on the large e-pockets in the d-wave state, 
due to the hybridization between two e-pockets that is 
inherent in 122 systems, which is inconsistent with the 
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FIG. 2: (a) Qs- (Qc-) dependence of Ae for d-wave (s++-wave) 
state at T = 0.03 eV. (b) SC gap functions for S-|-+-wave state. 



specific heat measurements^'^ that report the isotropic 
SC gap. We propose that the study of impurity effect 
on Tc is useful since nodal d-wave state is fragile against 
impurities. 



III. ORTHORHOMBIC STRUCTURE 

TRANSITION AND SOFTENING OF SHEAR 

MODULUS Cs 

A. Tviro-orbiton mechanism 

In many Fe-based superconductors, the SC state oc- 
curs next to the structure transition, suggesting a close 
relation between these two phenomena. Since the lat- 
tice deformation in the orthorhombic phase is very small 
((a — 6)/(a -!-&)< 0.003), we can safely rule out the 
cooperative Jahn- Teller effect due to electron-acoustic- 
phonon. Hereafter, we promise that a = xaxis and b — y 
axis. Below, we discuss the structure transition due to 
ferro quadrupole-quadrupole interaction. 

For this purpose, we introduce the strain-quadrupole 
coupling: 



Hs = Vs^ es<i 



(5) 



where es is the strain for the orthorhombic deformation 
(es oc a — 6), 0s is a quadrupole operator that belong to 
the same representation of eg, and 775 is the coupling con- 
stant. Fernandes et al.^'^ considered the "spin quadrupole 



operator" (j)s = si -52, where si represents the spin op- 
erator at sublattice I = 1,2. (In their theory, stripe- type 
commensurate magnetic correlation is assumed.) In this 
case, (ps is a non-local operator. On the other hand, the 
present authors studied the "charge quadrupole opera- 
tor" (j>s — 0^2 -y2, which is a local operator. In the pres- 
ence of the strain-quadrupole coupling, the shear modu- 
lus Cs is given as 



Cs'^Cs},a + gsXsiO,0)) 



(6) 



where Csfl is the lattice shear modulus, and gs — 
Vs^so- Xs{<i,^) is the total quadrupole susceptibil- 
ity, given by the Fourier transformation of xs i^^ t) = 

(T.0q(r)0_q(O)). 

Experimentally, Cs follows the Curie- Weiss relation in 
the under-doped systems with T5 > 0: 



Cs 
Cs,o 






(7) 



On the other hand, Cs deviates from the Curie- Weiss law 
in the over-doped systems without structure transition. 

From now on, we consider the charge quadrupole sus- 
ceptibility for (/) = Ox2_y2, x^2_„2(0), based on the or- 
bital fluctuation theory. In the RPA, the quadrupole in- 
teraction (1) together with the Coulomb interaction give 
rise to the AFQ fluctuations for the channels T = xz, yz 
and xy, while Xx^-y^iO) remains small. Therefore, Cs 
softening cannot be explained within the RPA. How- 
ever, the mode-mode coupling process with respect to 
Xxzi'i) gives the development of ferro-quadrupole (FQ) 
fluctuations; X%.y2i0) - T j:^{xgici)}\ This pro- 
cess is called the "two-orbiton term" in Ref.^^ since the 
condensation of composite bosons with zero momentum, 
Oxzi^)Oxzi~<i)i induces the tetragonal structure transi- 
tion, because of the relation 0x2 -y2 ex O^^ — Oy^. The 
two-orbiton term has close similarity to the Aslamazov- 
Larkin theory of excess conductivity given by supercon- 
ducting fluctuations. 

By taking the two-orbiton term, the irreducible sus- 
ceptibility is given as 



xL"-y4q) 



X 



(0) 



!(q) +XTo(q) 



(8) 



where xi2„ 



j(q) is the bare bubble term, and its 



temperature dependence is weak. Within the classi- 
cal approximation, the two-orbiton term is XTo(q) — 
TA^^[^Xx2(q + l^)x?z(k), where A is the three-point 
vertex'^^. The exact expression is given in eq. (61) in 
Ref.26: 
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(9) 
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A^ 



where ^l){x) is di-Gamma function, and X 

[c was introduced in eq. (2).) Since XTo(q) oc T^^ for 



An 



q — 7- 0, the two-orbiton term shows strong T-dependence. 
When [7 = 0, the total quadrupole susceptibihty is given 



x^2„,2(q) = xL"-,=(q)/(i - 5sxL"^,2(q)) (lo) 

where gs is the quadrupole-quadrupole interaction due 
to acoustic phonon. As a resuh, the elastic constant Cs 
is given by eqs. (6), (8) and (10), together with eq. (61) 
in Ref.26. 
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FIG. 3: Cs/Cs,a given by the two-orbiton theory that repro- 
duce experimental data for Ba(Fei_a;Coi)2As2 observed by 
Yoshizawa . 



theoretically, the obtained Tafq slightly exceeds Ts for 
X — 7 ^ 8%, because of the lack of experimental data 
at low temperatures. In fact, ^^ is expected to deviate 
from the Curie- Weiss behavior at low temperatures, as 
predicted by the SCR theory. In Fig. 4, is obtained 
experimentally by using the relation (1 — Cs/C's.q)^^ oc 
T — 6 when Cs/Cs o follows the Curie- Weiss relation in 
eq. (7). 
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In Fig. 3, we plot Cs/Cs, o based on the two- 
orbiton theory. We put orbital fluctuation parame- 
ters as I = 1200K, I' = 1400K, and Tafq = 92K- 
— 120K. These values are consistent with theoretical cal- 
culation based on the FLEX approximation^^. As for 
the spin fluctuation parameters in cuprates, I — 0.2eV 
(O.leV) for Lai_2;Sra;Cu02 (YaBasCusOy). We also 



X 



(0) 



.(0) 



2.88eV"\ and 



set gs = 0.183eV, as 

X = XTo{0)/eT^5A7^7.u: 

In the FLEX approximation^^, Tafq changes from pos- 
itive to negative by carrier doping, while other parame- 
ters (such as I and I') are insensitive to the doping. We 
can fit the recent experimental data by Yoshizawa et al.'^ 
for Ba(Fei_a;Coa;)2As2 with x = ~ 16%, by changing 
Tafq from 90K to -120K, together with X = 5.47 ^■ 
7.11. The fitting data shown in Fig. 3 reproduce the 
experimental data of Ref."* almost perfectly. This fact is 
a strong evidence for the success of orbital fluctuation 
theory in iron pnictide superconductors. 



B. Phase diagram and SC transition temperature 

We show the phase diagram for Ba(Fea;Coi_2,)2As2 in 
Fig. 4. Tc and Ts shows the experimental values, and 
Tafq is obtained by the fitting using the two-orbiton 
theory. While Tafq should almost coincide to Ts at Ts = 



FIG. 4: Phase diagram for Ba(Fea;Coi_2:)2As2 obtained by 
the orbital fluctuation theory. Tafq is derived from the tem- 
perature dependence of Cs observed in Ref.*. We also show 
Tc derived from eq. (12) by solid squares using orbital fluc- 
tuation parameters derived from Cs 

Based on the obtained orbital fluctuation parameters, 
we calculate the SC transition temperature. Here, we 
consider the orbital-fluctuation mediated s+_|_-wave state 
due to inter-pocket nesting. The linearized gap equation 
is given as 

AijAa(fc,e„) = TTT^y'V' / ^Va,/3(fcen,pem) 



j3,p m 



(11) 



where Xe is the eigenvalue; As = 1 is satisfied at T = Tc. 
a,/3 represent Fermi surfaces (FSs), k[p) is the momen- 
tum on the FS a (FS /?), and Va^p{ken,pem) is the 
interaction between (a, /c) and (/3,p). Here, we con- 
sider the attractive interaction due to orbital fluctuation; 
Va.j3 OC Xxz(vz)- ^^^ simplicity, we assume that both the 
Fermi velocity and SC gap are isotropic. Then, the mo- 
mentum integration along the FS in eq. (11) can be per- 
formed as /„ ^Va,i3{ken,pem) = C^/(l + |e„ - e„|/a;o). 



where C is a constant. If we assume A^ — A^, we obtain 
the following simplified gap equation: 



C. Difficulties in other orbital fluctuation theories 



^^A(.„) = r5:5:- 



Ci-ix* 



A(e™) 



fi m 



1 1/^0 



(12) 



where [i* is the Morel- Anderson pseudo potential; we put 
^x* = 0.1. We put C = 0.075 to obtain Tc = 25K when 
^AFQ = 0. In Fig. 4, we show the SC transition tem- 
peratures obtained by eq. (12) by solid squares, which 
are similar to experimental values. This fact means that 
the orbital fluctuation parameters derived from Cg are 
reasonable. 
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FIG. 5: Resistivity given by eq. (13), using the parameter 
Tapq derived from Cg. We can reproduce experimental data 
in Ba(Fei_a;Coa;)2As2^^ below 200K very well. 



Next, we discuss the non- Fermi-liquid behavior of the 
resistivity caused by the AF orbital fluctuations. In 
Fig. 5, circles show the resistivity of single crystal 
Ba(Fei_a;Coa;)2As2 reported in Ref.'^^. According to 
spin/orbital fluctuation theory^^, the resistivity in two 
dimension is given as 



AT-^e 



Po 



(13) 



where ^ is the correlation length of AF fluctuations, po 
is the residual resistivity, and ^4 is a constant. Using this 
equation together with ^^ = 1/{T — Tafq) derived from 
Cs, we can reproduce experimental data below 200K very 
well, as shown in Fig. 5. Here, we put A = 0.49 for 
all X, and (rAFQ,Po) = (-8,80) for x = 8%, (0,72) 
for X = 10%, (-50,70) for x = 12%, and (-90,67) 
for X = 14%, respectively. (Tafq's are obtained from 
Fig. 4.) In the present theory, Fe-ion optical phonons 
together with the Coulomb interaction induce the AF- 
orbital fluctuations, which give rise to the s++-wave state 
and non-Fcrmi-liquid-like transport phenomena. More- 
over, two AF-orbitons with zero total momentum in- 
duce the ferro-orbital fluctuations, which are the origin 
of the orthorhombic structure transition. Therefore, the 
present orbital fluctuation theory can explain both the 
structure transition and the superconductivity. 



Here, we discuss various difficulties in other or- 
bital fluctuation models proposed for iron-based 
superconductors^^''^'^ To explain the orthorhombic struc- 
ture transition, the ferro-Oa;2_y2 fluctuations should be 
the most divergent. The U' > U model in Ref."^^ cannot 
explain the structure transition since the most divergent 
fluctuation is 0^2 -type. 

Next, we consider the "Oj.2_j,2-fluctuation model, 
in which all quadrupoles other than 0x2 _y2 do not 
fluctuate. Then, electrons with xy-orbital character, 
which occupies one-third of the total DOS at Fermi 
level, are not involved in the 02,2_y2-fluctuations, since 
{xy\Ox2-y2\m) = for any d-orbital m. This fact means 
that "gapless SC state with large residual DOS" is re- 
alized in the 03,2 _y2 -fluctuation model. Figure 6 shows 
the gap structure given by the (a) "O^-^/jj^-fluctuation 
model" in Ref.^^ and (b) "Oj.2_y 2 -fluctuation model" in 
Ref.'^^. In the latter model, the quadrupole interaction 
H' = -9{'^)J2iOl2_y2 ■ Ol2_y2 is introduced, consid- 
ering the As-ion acoustic phonons. In (a), fully-gapped 
s++-wave state is realized. In (b), in contrast, the SC 
gap on the xy-orbital hole-pocket around (tt, tt) in the 
unfolded BZ is almost gapless, and the SC gap on e- 
pockets is highly anisotropic. Therefore, 0^2 _y2 fluctua- 
tion model cannot explain the fully-gapped nor nodal gap 
structure observed in almost all iron-based superconduc- 
tors. 

We also comment on the "T-linear resistivity" observed 
in various iron-based superconductors. In the ferro- 
Oj.2_y2 fluctuation model, the portion of Fermi surfaces 
with xy orbital would become the cold-spot with 7 ex T^. 
Thus, a conventional Fermi liquid behavior p oa T^ would 
be obtained in the ferro-Oj.2_y2 fluctuation model, incon- 
sistently with experiments. 

Therefore, orbital fluctuations theories in Refs.'^^''^^ 
have serious difficulty in explaining the "fully-gapped" 
s_|__|_-wave state as well as non-Fermi liquid transport phe- 
nomena in iron-based superconductors. 



D. Eff'ective Action 

In Ref.^^, we have calculated the shear modulus Cg 
using the Green function method. Since the derivation 
was rather complex, we rederive the same expression 
in a more simple manner by introducing the Hubbard- 
Stratonovich field, in analogy to the analysis in Ref.'^'*. 
By performing the Hubbard-Stratonovich transforma- 
tion, we introduce the following effective action S = 
Sq + Sint + •S'stiain for the charge quadrupole field (fir 
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FIG. 6: The SC gap functions for a = 0.98 in the (a) O^z/yz- 
fluctuation model and (b) 0^2_y2-fluctuation model. In (b), 
the a::2/-orbital hole-pocket around (tt, tt) is almost gapless, 
since a;y-orbital is not involved in the 02,2_y2-fluctuations. 
In (a), the eigenvalue of the gap equation (11), A^, is larger 
with smaller interaction g. 



X^2_y2{q) 
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= A^ 
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(18) 



where / • • • — TJ2n I JW^ ' " ■ Here, we have dropped 
both the vertex and self-energy corrections. Apart from 
the factor A^ , the first term the irreducible susceptibility 
(18) is equivalent to xto in eq. (9), or xo,ncm in Ref.'^^. 
We also added the bare susceptibility y^ 2_ 2 as the sec- 
ond term in eq. (18). Due to the two-orbiton term, 
X"^"^ (0) is strongly enhanced near AFQ-QCP in propor- 
tion to ^. When 7;^ is finite, x 2_ 2(0) diverges even 
if ^ is finite, and therefore the orthorhombic structure 
transition takes place in iron-based superconductors. 

Note that the present action in eqs. (14)-(16) are 
mathematically equivalent to eqs. (3) and (4) in the 
spin-quadrupole theory in Ref.^^, by replacing 0r (r = 
XZ.YZ) with 0i (i = 1,2), and Xr('?) with x'^A)- 
Thus, the derived T-dependences of Cs are essentially 
the same. However, the coupling between strain and 
spin-quadrupole, 775, would be too small to fit experi- 
mental data: We will discuss this issue below. 



{V = XZ,YZ): 
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] = ^ E /{x?('?)}"V^(g), (14) 



-f {4>xz{x)4>Yz{x)Y (15) 
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J X 

where J^' ' ' = Jo dr J cP x ■ ■ ■ , So represents the most 
divergent AF quadrupole susceptibility Xxz(yz)^1^^ "^int 
is the quadrupole-quadrupole interaction due to acoustic 
phonons, and S'strain is the strain-quadrupole coupling. 
In this subsection, we introduce the XF-coordinate that 
is —45 degree rotated from the a;y-coordinate along z- 
axis. Apparently, (l)xz = {4'xz ~ (l)yz)/V^, (f^YZ = {4>xz + 
(f)yz)/^, and Ox-2_y2 — OxY- In Eqs. (15) and (16), 
g'g = gs ■ A^ and rj'g = rjs ■ A. where A is the three-point 
vertex with respect to {Oxy,Oxz,Oyz) that had been 
analyzed in Ref.^^. 

Now, the quadrupole susceptibility X^2_,2 is given by 

X y 

the second derivative of the partition function with re- 
spect to r]es{x). If we put g'g = 0, the quadrupole sus- 
ceptibility is given as A'^Xxz(^)XYzi~^)- ^y performing 
the Fourier transformation and taking the interaction g'g 
into account, the total quadrupole susceptibility is ob- 



E. Comparison with spin quadrupole theory 

As discussed in Sec. Ill, the orthorhombic structure 
transition is cause by the divergence of the FQ suscepti- 
bility with i^ — y'^ symmetry. In this paper, we discussed 
the FQ fiuctuation induced by the AF quadrupole fluctu- 
ations, due to the two-orbiton process. That is, antiferro- 
orbital fiuctuations with respect to Oxz £^nd Oyz (=or- 
bitons) induce the s++ wave superconducting state, and 
the bound-state formation of two-orbitons with zero mo- 
mentum, Oxz(q)Oi'z(— q) ^ Oxy(O), give rise to the 
development of x$y(0) = X^2_y2(0). 

In this subsection, we discuss the "spin-nematic the- 
ory" for the structure transition discussed in Refs.'^"'''^*, 
and explain a close relation to the two-orbiton mecha- 
nism. Their theories can be interpreted as the "two- 
magnon process" in our language. They had studied the 
non-local spin quadrupole operator cjjs = ini ■ m2 as ex- 
plained in Sec. Ill A, and shown that the Ising-like order 
0s 7^ In both theories, the Ising-like order ((0s) 7^ 0) 
occurs prior to the vector order {{nii), {{Oxz, Oyz)) ^ 0), 
since the latter is easily suppressed by thermal and quan- 
tum fluctuations. When the boson is orbiton (magnon), 
the realized superconductivity is the s++ (±) wave state. 

According to the fitting done in Sec. Ill, we have to 
assume gs ~ O(0.1)eV to reproduce Ejt = Ts — ^ 
O(10)K; the corresponding dimensionless e-ph coupling is 



A < 0.1. In this case, the required strain-quadrupole in- 
teraction 775 is about 0.5eV/ Angstrom. As for the strain- 
charge-quadrupole interaction, 77^ ^"^^^ ^ 0.5eV is actu- 
ally obtained by the point-charge model^^. On the other 
hand, the strain-spin-quadrupole interaction ry^"" would 
be of order ~ 6J{R)\m\'^/dR, where J{R) is the nearest- 
neighbor magnetic interaction and \m\ is the magnetic 
moment. Since 5J{R)/5R ^ 0.04eV/ Angstrom accord- 
ing to the first principle study"^^, we expect that 775"" 
is one order of magnitude smaller than 77^ ^'^^'^. Then, 
there seems to be difficulty in reproducing the softening 
of Cs that occurs for a wide temperature range (from 
the room temperature to Tg; see Fig. 3) in terms of the 
spin-nematic scenario. 

Finally, we note that the spin-quadrupole scenario re- 
quires that Xs{^) is "commensurate". However, recent 
neutron experiment revealed that the magnetic order be- 
low Ttv (< Tg) in Ba(Fei_a;Co2;)2As2 is incommensurate 
for X > 5.7%, although structure transition occurs for 
X > 7%. This fact would support the realization of two- 
orbiton mechanism, since it does work even if Xxz( z')(*^) 
is incommensurate. 



IV. IMPURITY-INDUCED NEMATIC ORDER 
A. mean-field approximation 

In previous sections, we have shown that both s++ 
wave superconductivity as well as the structure transition 
originate from the AF orbital fluctuations, xiz(vz)^Q) 
caused by the inter-orbital nesting. In this section, we 
discuss the nematic ordered state in the tetragonal phase, 
which had been observed in many pnictides. For ex- 
ample, in Ba(Fei_2;Co2;)2As2, large in-plane anisotropy 
starts to occur about lO^lOOK higher than Ts^ ■ This 
"nematic electronic state" free from lattice deformation 
had been also observed by ARPES measurements^'^''*'^'^^, 
optical conductivity measurements''^, and the anisotropy 
in the magnetic susceptibility"^. The origin of this ne- 
matic order state had been one of the great open issues 
in iron-based superconductors. 

Recently, local-density-of-states (LDOS) around im- 
purity sites had been studied in detail by STM/STS 
technique, and nontrivial breakdown of the four-fold 
symmetry (C4^) had been found in many pnictides. 
For example, single Co impurity induces a remark- 
able non-local change in the DOS along the a-axis 
(longer axis) with a length of ^ Sapo-Fo in under-doped 
Ba(Fei_a;Coa;)2As2^"^. Such a drastic impurity-induced 
anisotropic state is nontrivial since the orbital polariza- 
tion in the orthorhombic phase is very small {rixz — nyz ~ 
O(10~^)). Similar impurity-induced non-local change in 
the DOS was also observed in LiFeAs with tetragonal 
structure by Hanaguri et at^^. In many cases, impurity- 
induced states break the C4t,-symmetry, belonging to 
C'2v, (^2, or Cih- 



It would be natural to expect that impurity-induced 
low symmetric state explains the nematic phase above 
Ts- Recently, in-plane resistivity had been mea- 
sured in detwinned Ba(Fei_j,Co2:)2As2 crystals with high 
quality '*^''*^ . They found that pb — pa is proportional to 
the Co concentration a; for a; < 4%. Moreover, the resid- 
ual resistivity per 1% Co impurity reaches ~ lOO/iilcm, 
which exceeds the unitary limit of the local impurity po- 
tential (^ 2QpVlcTa). Therefore, it had been proved ex- 
perimentally that a Co impurity works as a strong scatter 
in the FeAs plane, like Zn impurity in cuprate supercon- 
ductors. 

In strongly correlated electron systems, impurity po- 
tential frequently causes nontrivial non-local change in 
the electronic state. For example, in nearly antiferro- 
magnetic metals, antiferromagnetic correlation is drasti- 
cally enhanced near the nonmagnetic impurity site, and 
the local magnetic moment emerges around the impurity. 
Such phenomena are indeed observed in under-doped 
cuprates^^. As for the iron-based superconductors, the 
system would be close to antiferro-orbital QCP. There- 
fore, it is natural to expect the occurrence of "impurity- 
induced non-local orbital order" in pnictides. 

In this section, we study the single impurity prob- 
lem with potential / in ten-orbital tight-binding model 
for pnictides, using the mean-field approximation in 
real space Since we concentrate on the impurity-induced 
orbital-ordered state, we take only the quadrupole- 
quadrupole interaction in eq. (1) into account: We con- 
sider that the coupling constant g in eq. (1) originates 
from both the e-ph interaction as well as the multiorbital 
Coulomb interaction, as discussed in Sec. II. We ob- 
tain the following mean- field self-consistently for g < gc, 
where gc ~ 0.221 is the critical value for the bulk orbital 
order"*^: 



^l,m — \Ci^iaCi,ma)l,g {Ciic,Ci,ma) 1 ,0 



(19) 



where i is the Fe site, and l,m represent the d-orbital. 
Note that M/^ is impurity-induced mean- field for g < 
gc since it vanishes when / = 0. Then, the mean-field 
potential is given as 



Z^ 
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(20) 



and the mean-field Hamiltonian is -ffMF — Hq + ^^ 5* + 
const. F"^ is the bare four-point vertex for the charge 
sector. We solve eqs. (19)-(20) self-consistently. 

Since the present mean-field has 15 components at each 
site, it is convenient to consider the following quadrupole 
order: 



01 



22^0r 



'M, 



l.m 



(21) 



Ot^ is the 



where F = xz, yz, xy, z , and x — y . uj. 
matrix element of the charge quadrupole operator given 
in Rcf.^*. We note that the hexadecapole {I = 4) order 
is negligible in the present study. 



In Fig. 7, we show the obtained DOS in real space for 
the cluster of 800 Fe sites with a single impurity site^^. 
First, we consider the case / = — 2eV in (a)-(c). For 
g = 0.200 (a), the impurity induced mean-field is absent. 
The small modulation of the LDOS around the impurity 
is caused by the Friedel oscillation. For g > 0.206, im- 
purity induced non-local orbital order is induced. The 
symmetry of the corresponding LDOS is C2v, as shown 
in (b) for g = 0.210. The suppression of the DOS is 
caused by the non-local orbital order, consistently with a 
recent optical conductivity measurement^^. The LDOS 
are further suppressed with increasing g, and the sym- 
metry is lowered to C2 for g > 0.212. In (c), we show 
the numerical result for g = 0.216 with C2-symmetry. 
Similar impurity-induced LDOS with C2(y)-symmetry is 
also obtained for a positive impurity potential. In con- 
trast, we obtain the orbital order with C4 symmetry for 
/ = -t-leV. 

(a) l=-2eV, g=0.200eV (b) l=-2eV, g=0.208eV 
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FIG. 7: Obtained LDOS given by the mean-field theory in 
the presence of impurity potential at (0,0). (a) (/,<?) ~ 
(-2,0.200): without orbital order, (b) {I,g) = (-2,0.208): 
orbital order with C2v-symmetry. (c) (/,(/) = (—2,0.216): 
orbital order with C2-symmetry. (d) (/,<?) = (-1-1,0.216): or- 
bital order with C4-symmetry. 



The dominant impurity-induced quadrupole orders are 
antiferro Oxz^ Oyz and O^y, reflecting the quadrupole- 
quadrupole interaction in eq. (1). The indispensable 
ingredient for C2-order is the quadrupole interaction 
for F = xz/yz channels. We have verified that the 
quadrupole interaction for F = a;^ — y^ channel, which is 
caused by acoustic phonon, cannot realize the C2-order. 

By symmetry, the obtained C2-order can be aligned by 
the strain-induced quadrupole potential, which is given 
in eq. (5) with ^5 = Ox2-y2. The effective quadrupole 



potential is A^ = rises ■ X^2_j,2(0)/x^°2_j^2(0), which is 
strongly enhanced near Ts due to the two-orbiton pro- 
cess as discussed in Sec. Ill A. This would be the rea- 
son why the nematic ordered state is easily detwinned 
by small uniaxial pressure near Ts- In fact, detwinning 
by uniaxial pressure is possible only when the structure 
transition is the second-order^^, in which case the factor 

X%-y2iQ)/x°2_y2iO) develops divergently. In over-doped 
systems, detwinning by pressure would be difficult since 
x3^_^2(0)/xi°2^_ 2(0) is no more large^^. 



B. in-plane anisotropy of resistivity 

According to ARPES measurements in detwinned sys- 
tems, AE < for a; = a-axis is longer than y = b- 
axis, i.e., Uxz > nyz^^''^^. Then, experimentally observed 
(tt, 0) SDW order should be induced theoretically^^. In 
this case, the C2 order aligns along the x = a-axis in 
the present numerical results in Figs. 7 (c) and (d). We 
have calculated the in-plane resistivity pa and pb in the 
presence of dilute C2 orders along a-axis using the T- 
matrix approximation, which gives the exact result when 
the impurity concentration is dilute and localization is 
negligible. We find that pa is smaller than pt, by ^ 40%, 
consistently with experimental reports. In Fig. 8 (a), 
we show the alignment of impurity-induced C2 objects 
under uniaxial pressure. The anisotropy of resistivity in 
the nematic phase given by the T-matrix approximation 
is shown in Fig. 8 (b). 

Here, we explain how to calculate the resistivity in 
the presence of aligned non-local orbital orders (nematic 
phase). The T-matrix is given by solving the following 
equation in the orbital-diagonal basis: 

Tr.r'{oj) — (I + S)rSr,r' 

+ J2ii + S)MX"i^)T^'yi^) (22) 



where S is the mean-field potential, and Gr (w) is the 
Green function in real space without impurities. Ir = 
IlSr,o is the impurity potential. Note that Sr is diag- 
onal with respect to the position since our interaction 
Hamiltonian is local. 

The preset T-matrix is non-local because of Sr for r 7^ 
0. After the Fourier transformation, the self-energy due 
to multiple scattering in the T-matrix approximation is 
E(k, oj) = nimpTk,k(w), and the full Green function is 
G{k,Lo) = {lu + p- H°~ t^{k,uj))-^. Then, the in- 
plane conductivity is given by 



^E 



v^^JUGciK^SW 



(23) 



where v = x,y, and a represents the ath band. Uk,i/ is 
the group velocity of the ath band, and Gaik^uj) is the 
full Green function in the band-diagonal basis. Jk.,y is 
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the total current including the current vertex correction, 
which is given by solving the following Bethe-Salpeter 
equation: 



ja 



where Ik k' 



^k,iy 



E4"ilGMp,^^)P^p^ 



(24) 
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As shown in Fig. 8 (b), we obtain pt/pa ~ 2 in the 
nematic phase: The anisotropy is enhanced by including 
the current vertex correction. 

Note that the averaged residual resistivity {pa + Pb)/'^ 
per 1% impurity with / — +leV or — 2eV reaches 
^ 50/ir2cni, which is comparable to the residual resis- 
tivity by 1% Co impurities in Ba(Fei_a;Coa;)2As2; 50 ^ 
100/if}cm".46 



(b) 
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C. Comparison between theory and experiments 

In a pure two-dimension system, a single nonlocal or- 
bital order in the mean-field approximation might disap- 
pear by thermal and quantum fluctuations. In real sys- 
tems, however, orbital order would be stabilized by the 
|T^j^,(i(5)p is the irreducible vertex. correlation between impurities for riimp > 1%- Since the 

AF fluctuations increase as T decreases in real systems, 
we can interpret that g monotonically increases as T de- 
creases in the present mean-field approximation. There- 
fore, we expect that impurity-induced non-local orbital 
order is stabilized below the nematic transition temper- 
ature T*, at which g{T) would be close to gc- 

Although orbital order is absent above T*, strong or- 
bital fluctuations should appear near the impurity site, in 
analogy to the impurity-induced strong magnetic fluctu- 
ations realized in under-doped cuprates^^. Therefore, in- 
elastic scattering given by impurity-induced strong fluc- 
tuations causes large "residual resistivity" even above the 
impurity- induced Neel temperature^^. By the same rea- 
son, residual resistivity in under-doped pnictides would 
be large even above T*. The order parameters for the 
impurity induced nematic order are antiferro-0^2/j^2,/2,j^, 
while the order parameter for the structure transition is 
fcrro-Oa:2_j,2. Because of the difference in order param- 
eters, Ts would be rather insensitive to the presence of 
impurities. 

The nematic transition also occurs in the tetragonal 
phase of BaFe2(Asi_a:Pa;)2, which was confirmed by the 
in-plane anisotropy in the magnetization [xa 7^ Xb) us- 
ing the torque measurement under magnetic field^. In 
BaFe2(Asi_a;Pa;)2, P sites would works as impurities, 
which give finite potential on the neighboring for Fe sites. 
In this case, we had verified that non-local orbital or- 
der with C2- or Ci/i-symmetry appears around the im- 
purity site in the mean-field approximation. In contrast, 
nematic order had not be observed in (K,Ba)Fe2As2^^, 
maybe because impurities outside of FeAs plane would 
be too weak to induce orbital order. 

In this paper, we consider that the charge quadrupole 
order occurs at Ts, and the nematic order at T* (> Ts) 
originates from the impurity-induced orbital order. How- 
ever, difference scenario had been proposed by the au- 
thors in Ref.^: They consider that the quadrupole order 
occurs at T*, and T5 is just a meta-transition without 
symmetry breaking. In their phenomenological Landau 
model, very small orthorhombicity (a — b)/{a + b) oc- 
curs below T* , and it increases drastically below Ts as a 
meta-transition, consistently with experiments. 
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FIG. 8: (a) Detwinning of the impurity-induced 

C2-orders under uniaxial pressure. The strain- 

induced effective quadrupole potential on each Fe is 
Vs(^s{x%_y2{0)/x^°2_y2(0))Ox2-y2. (b) lu-plaue anisotropy 
of resistivity in the nematic phase obtained by the T-matrix 
approximation. We assume 1% impurity with / = — 2eV. 

Here, we have studied the resistivity due to elastic scat- 
tering caused by spatially extended effective impurity 
potentials. In addition, strong antiferro fluctuations in 
pnictides should present various non-Fermi liquid trans- 
port phenomena in the normal state. In fact, relations 
p (X T and i?H oc T~^ were observed in BaFe2(As,P)2^°. 
Although such behaviors are frequently ascribed to the 
evidence of spin fluctuations^-'^, they are also brought by 
the development of antiferro-orbital fluctuations^®. 



V. DISCUSSION 

In the present paper, we have studied the realistic five- 
orbital model for iron pnictides. It was found that the 
Oajz-AFQ fluctuations develop owing to the quadrupole- 
quadrupole interaction due to e-ph and Coulomb inter- 
actions. This fluctuations not only cause the s-|__(--wave 
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superconductivity, but also the orthorhombic structure 
transition due to two-orbiton process. Using the two- 
orbiton term, we can fit the recent experimental data 
of Cs in Ba(Fei_2;Coa;)2As2^ for wide range of doping, 
only by choosing Tafq while other parameters are almost 
fixed. This fact is a strong evidence for the success of or- 
bital fluctuation theory in iron pnictide superconductors. 
However, we consider the impurity effect on the s_|_+- 
wave state: In the weakly correlated metals, impurity 
effect on the s_|__|_-wave state is very small, known as 
the Anderson theorem. In strongly correlated metals, 
in contrast, impurity-induced change in the many-body 
electronic states could violate the Anderson theorem. 
In Sec. IV, we have shown that the impurity-induced 
orbital order causes remarkable reduction in the DOS, 
accompanied by the suppression of orbital fluctuations. 
For this reason, orbital-fluctuation-induced s++ wave 
state should be suppressed by impurities in the nematic 
state. In Ba(Fei_a;Co2:)2As2, the suppression in T^ per 
1% Zn-impurity, which gives a strong impurity poten- 
tial, is -ATc/% - 3Ki^ while -AT^/% ~ 20K is ex- 
pected in the s±-wave state when the mass-enhancement 
is m*/mb ~ 3^'^. Such small impurity-induced suppres- 
sion of Tc in Ref.^^ would be naturally explained by the 
impurity-induced orbital order. 
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FIG. 9: The phase-diagram for iron-pnictide superconductors 
obtained by the present orbital fluctuation theory. Ts is the 
orthorhombic transition temperature (= FQ order tempera- 
ture), and Tn is SDW transition temperature. The fact that 
two QCPs at Ts = and Tafq — almost coincide means 
that novel "multi orbital QCPs" are realized in iron pnictides. 
At Tafq, the AFQ-order does not occur since it is prevented 
by the FQ order at Ts. The impurity-induced nematic state 
is realized below T*. 



In Fig. 9, we summarize the phase-diagram of iron- 
pnictides given by the present orbital fluctuation theory. 
We stress that Tafq, which is determined experimentally 
from Cs, is positive in the under-doped case {Ts > 0) 
while it is negative in the over-doped case. Especially, 
Ts ^ for Tafq = 0, consistently with experiments^'^. 
This result indicates that QCPs for AFQ and FQ or- 
ders almost coincide at the endpoint of the orthorhombic 
phase. Therefore, AFQ-QCP is not hidden inside of the 
orthorhombic phase, favorably to the orbital-fluctuation- 
mcdiatcd s+_|_-wave SC state^^"^^. Below T*, impurity- 
induced orbital order with C2(\/)-symmetry starts to oc- 
cur. Each orbital order can be aligned by applying 
tiny uniaxial pressure, resulting in the large in-plane 
anisotropy in the resistivity. 



In summary, the present orbital fluctuation theory can 
explain the (i) superconductivity, (ii) orthorhombic tran- 
sition accompanied by large softening of Cs, and (iii) 
impurity-induced nematic order. To explain (i)-(iii), 
we have to assume that Xxz(vz) ^^ ^^^ most divergent 
quadrupole susceptibility, which is actually satified in 
the present model Hamiltonian. We have derived the 
orbital fluctuation parameters from Cs, and succeeded 
in explaining Tc and p(T) using the derived parameters. 
These results are strong evidence for the realization of 
the orbital-fluctuation-mediated s_|--)_-wave superconduc- 
tivity in iron pnictides. 
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